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SUMMARY 


A method Is presented for studying the transient "behavior of the 
flapping motion, as well as for calculating the steady-state flapping 
an 5 )lltudes , of free-to-cone and seesaw rotors operating at extreme -Plight 
conditions. The method is general and can "be applied to "blades of any 
airfoil section, mass dlstri"btctlon, twist, plan-form taper, root cutovrt, 
and flapping-hinge geometry. Stall and compressibility effects can also 
be accounted for. 

The method is Illustrated by exangoles of blade-stability calcula-. 
tions of free-to-c.one and seesaw rotors at tip-speed ratios equal to or 
greater than 1.0 and of the transient flapping response -of a free-to-cone 
rotor to arbitrary control irqiuts at a tip-speed ratio equal to 0.3. 


INTRaDUCTIQN 


The flapping motion of autogiro and helicopter rotor blades has 
shown Itself, both in theory and actual practice, to be -very stable for 
conventional tip-speed ratios (that is, below about 0.5). Some doubt 
exists, however, as to the stability of blade motion at tip-speed ratios 
equal to or greater than 1.0, a question that is of interest in connec- 
tion with the "unloEided rotor" type of helicopter operation. At these 
extreme tip-speed ratios,- the aerodynamic and blade -inflexibility assump- 
tions en^jloyed in existing analyaes of the problem (in refs. 1 to 5, for 
exanple) became questionable. The aerodynamic factors introducing the 
most uncertainties in these analyses are the assumptions of -\installed 
blade sections and Inconpresslble flow, the neglect of the reversed- 
veloclty region, and the use of small— angle assumptions in connection with 
the section inflow and blade flapping angles. 
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It is the object of this paper to develop equations and to present 
a method for determining the stability of blade flapping motion vhlch 
avoids the no-stall and incompressibility assumptions, as well as some 
other restrictive aerodynamic assumptions that are normally employed in 
rotor analyses, and which properly accounts for the reversed-velocity 
region. The method, which is applicable to numerlc&Q. solution by automatic 
computing machines, can also be iised to account for the effects of stall, 
compressibility, emd high rotor angle-of -attack operation on the steady- 
state flapping motion at more conventional helicopter tip-speed ratios 
and to compute the transient blade-flapping response arising from arbitrary 
control inp\its. No specific terms are included in the equations, however, 
to account for blade torsional or bending flexibility. The method is 
illustrated by examples of blade-stability calculations of free-to-cone 
and seesaw rotors at tip-speed ratios equal to or greater than 1.0 and of 
the transient flapping response of a free-to-cone rotor to arbitrary con- 
trol inputs at a tip-speed ratio equal to 0.3« 


SYMBOLS 


A^,Bi 


B 




dD 

dPc 


coefficients .of -cos \lr and -sin ilr, respectively, in 
expression for 0; therefore, lateral and longitudinal 
cyclic-pitch angles, respectively, deg 

tip-loss factor (assimned equal to 0.97 herein); blade ele- 
ments outboard 'Of radius BR are assumed to have profile 
drag but no lift 

blade section chord, ft 

section profile drag coefficient 

equivalent blade chord (on thrust-moment basis), 

/ BR /r®® 

CT^drlJ r^dr, ft 

^c / 

section lift coefficient 

drag force on blade element, lb 

centrifugal force on blade element, lb 


dFj inertia force on blade element, lb 

dL lift force on blade element, lb 
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thruBt on blade element ^ lb 


offset of center line of flapping hinge from center line of 
rotor shaft j ft 

gravitational acceleration, ft/sec^ 

mass moment of inertia of blade about flapping hinge. 


/: 


[r - e) dr, slug-ft 


mass of blade per foot of radius, slugs/ft 

centrifugal-force moment of blade about flapping hinge, 

Ib-ft 

inertia-force moment of blade about flapping hinge, Ib-ft 
thrust moment of blade about flapping hinge, Ib-ft 

vel^t moment of blade about flapping hinge, Ib-ft 

wei^t moment of blade about flapping hinge at zero flapping, 

pR 

/ mg(r - e)dr, Ib-ft 
e 

helicopter rolling velocity, radians /sec 
helicopter pitching velocity, radians /sec 

distance measured along blade from axis of rotation to blade 
element, ft 

blade radius, ft 

time, sec 

resultant velocity perpendicular to blade-span axis at blade 
element, ft /sec 

conrponent at blede element of res\altant velocity perpendicu- 
lar both to blade-span axis and IIj , ft /sec 

s 

component at blade element of resultant velocity perpendicu- 
lar to blade-span axis and to shaft axis, ft /sec 
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induced inflow velocity at rotor (positive downward), ft /sec 

velocity along f 1 i ght path, ft /sec 

ratio of blade-element radius to arotor-blade radius, r/R 

angle between shaft axis and plane perpendicular to fli^t 
path, positive 'vdien axis is pointing rearward, deg 

blade-element angle of attack, meas\n:%d from line of zero 
lift, deg 

blade flapping angle with respect to sheift at particular 
azimuth position, radians 

first and second derivatives of ^ with respect to time t 

first and second derivatives of p with respect to azimuth 
angle t|t 

mass constant of rotor blade, CgpR^^Iji 

collective-pitch angle at blade root, average value of 
Instantaneous blade-root pitch angle around azimuth, deg 

collective-pitch angle at 0.75 blade radius, deg 


difference between root and tip pitch angles, positive when 
tip angle is larger, deg 

Instantaneous blade-section pitch angle; angle between line 
of zero lift of blade section and plane peipendicular to 
rotor shaft, cos sin ijr, deg 

V sin a- - V 

inflow ratio, 

flR 

V cos ttg 
OR 


tip-speed ratio. 



MCA TN 3566 


5 


p mass density of air, slugs/cu ft 

^ inflow angle at tlade element in plane perpendicvilar to 

"blade-span axis, tan"l ^®S 

i|r "blade azimuth angle measured from downwind position in 

direction of rotation, deg 

t 

SI rotor angular velocity, radians /sec 

Su"bscrlpt : 

c radius of cutout, that is, radius at which lifting surface 

of blade begins 


METHOD OF AMLTSIS 
Reference-Axis System 


In the analysis for the free-to-cone rotor as developed herein, the 
flapping hinge may he offset radially from the center line of the shaft. 
For such a configuration, flapping and feathering motions are not directly 
equivalent, and all angles and velocities are therefore referred to an 
axis coinciding with the center line of the rotor shaft. Also, even for 
zero hinge offset, the shaft axis is a convenient reference axis for cases 
in which the effect of arbitrary cyclic- or collective-pitch input on the 
blade transient motion is being stvidied. The rotor shaft is therefore 
used as the reference axis in this paper, although it sometimes may be 
more convenient to use the no-feathering axis as the reference, inasmuch 
as it is Independent of fuselage moment characteristics and center-of- 
gravlty location. 


Derivation of Blade-Flapping Equations 

The differential equation of blade flapping motion is derived for a 
free-to-cone rotor by equating the blade aerodynamic and mass moments 
about the flapping hinge. The origin of the various moments is shown in 
figure l(a), and the resulting eqmtlon is 


- Mg - Mj - M^j = 0 


( 1 ) 
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Expressions for the maments will now he derived. The centrifugal- 
force momenb is 


hut 


r® o r 1 

mfl^(r - e) e + (r - e)cos p sin p dr 


■ r® 2 

/ m(r - e) dr = Ij. 


( 2 ) 


(3) 


and 


r 


mg(r - e)cos p dr = = 1%'cos p 




(it should he noted that^ because Myf is sitiall in relation to the other 
hlade moments^ the cosine effect of a shaft tilt on has heen 
neglected.) After equations ( 3 ) and (4) are substituted into equation (2), 
the expression for centrifugal-force moment becomes 


Me » fl^sln p(% cos p + I 


(5) 


The egression for hlade inertia moment is 


nR 


Mj = / m|(r - e) dr 


- K (6) 

The Independent variable can he changed from t to i|r as follows: 


^ = ft M = np • 

dilf 

p = ft = ft^ i-| = ft^f 


dt dt 


dijr'^ 


( 7 ) 
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and equation ( 6 ) 130000168 




The expression for thrust momenb may be written with the aid of 
figure 1 as 

pBR 

Hji =J ipU%(r - e)cj cos 0 dr + 


1 2 

j ipU c(r - e)cfl^ sin 0 dr 


The profile-drag contribution to the thrust becomes significant primarily 
in the reversed-velocity region, at fli^t conditions involving high inflow- 
angles and hi g h tip-speed ratios. It is assumed^ of course^ that the cut- 
out r^j is equal to or larger than the' offse-t e, so that no moment is 

confuted for that part of the blade inboard of the flapping hinge. Let 


X = £ 

R 


I = e 
R 


Substituting equations (lO) into eqiuation (9) yields 


% = ipcfl^^ I u^(x - |)cj cos 0 dx + 


/ u^(x - sin 0 dx 
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The section nondimensional resultant velocity u and the inflow 
angle ^ may be found from the following equations: 


u = 



( 12 ) 


vtp^ = Ag cos p - (x - S)p - 


Hg sin p cos ilf + 


SL 

a 


X cos 


♦^5 


X 


sin '4f 


(15) 


= I + (x - I) COS P + Hg sin t 




(f) = tan“^ ( 15 ) 

(Note that, for the sake of sin 5 >licity, for terms involving pitching and 
rolling velocities in eq. ( 13 ) the helicopter is assxjmed to pitch and 
roll about the rotor hub.) 

The section lift coefficient c^ and drag coefficient c^.^ can be 

found from appropriate airfoil data. It should be noted that the method 
does not require an analytical representation of the section lift and drag 
characteristics. Thus, the airfoil characteristics corresponding to the 
particular blade section apd surface condition can be used, and stall, 

Mach nuniber, and Reynolds number effects can be accounted jfor whenever 
they" are deemed significant. 

T hfi section angle of attack can be found from the following relation: 


= 0 + ^ 


= 00 + cos sin '<1^+0 


(16) 
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In this eqiiatlon all the angles are in degrees. 

Substitirfcing equations (^),.(5)^ (8), and (ll) into equation (l) and 
dividing throu^ by yields 


1 pcR^ 

2 Iv, 


Z*® 2 /il.O p 

/ u (x - S)cj cos 0 dx + / u'^Cx - l)c^ sin 0 dx 


/ dMu* \ = Mrj'cos p 

sin p cos p + - p - -3L_ P = 0 

\ ^h / q2Ij^ 


By definition. 


r' 


pcR^ 


( 1 ?) 


(18) 


and equation (17) becomes 


|r- 


£ 


u2(x - |)cj cos JZJ dx + 


pl.O 

y^c 


2 

u (x - l)c^ sin ^ dx 


sin p fcos p + _ p _ Mw*cos p ^ ^ 

I " giby " n 2 


(19) 


If P is considered to be a small angle (that is, less than 20°), 
as it usual l y is, the equations can be same\diat simplified in that cos p 
can be assumed to be equal, to 1.0. It will also be noted that, althou^ 
a uniform-chord blade was assumed in the developnent of the thinist-moanent 
equation, in numerical work blades of any plan fom can be handled jiist as 
slmqoly by allowing c to remain beneath the integral sign. The blade mass 
constant 7', which is of value for blade-motion comparisons between 
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blades of different design, may then be retained by- basing 7' on an 

pm 

equivalent chord Ce equal to / 

Equation (19) is the basic differential equation describing the 
flapping motion of a free-to-cone rotor. As written, the aerod ynami c 
part of the equation is very general in that it can handle blades of any 
airfoil section, mass distribution, twist, plan-form taper, root CTibout, 

«Ti<^ flapping-hinge offset and C£in account for the effects of helicopter 
pitching or rolling velocities. The effect of flapping-hinge geometry 
where blade flapping or lagging motions result in pitch-a ng le changes can 
be ai-mply considered by the addition of the appropriate terms in equa- 
tion (16) . Stall anfl conqjressibility effects can be accounted for to the 
extent of actiial airfoil section data that correspond to the entire 

range of «w gl <=» of attack and Mach number encountered at all points of the 
disk, the reversed-velocity region. (The section Mach number 

is simply confuted as ufiR divided by the speed of sound.) No smaJ l .- 
ftn glft limitations have been made -regarding the blade inflow and flapping 
angles, althou gh it is assumed in the analysis that the blades are rigid 
in torsion and in bending and that radial flow effects can be neglected. 

The assun^tlon is also made that the rotational speed is a constant . Such 
an assxm^tlon is considered to be reasonable in that some form of rotor- 
speed control will probably be incorporated in most unloaded-rotor configu- 
rations because of the extreme sensitivity of rotor speed to smgle-of -attack 
changes at those conditions. 



The differential equation of blade motion for the seesaw rotor may 
be derived in an analogous manner, the only difference being that the 
moments of each blade aboiit the seesaw hinge are equated'. Thi» resulting 
equation is ' 



pi.o 

cos 0 dx -I- / u^c^Q sin 0 dx 

d^c 



cos 0 



sin 0 dx 




( 


2 sin p cos p - 



( 20 ) 



HACA TN 3366 


11 


where ijr and \jr + jt Indicate the Instantaneous azimuth positions of the 
two blades. It should he remembered that, in ^valuatlng equation (I3) 
for the seesaw rotor, the signs of ^ and 3 for each blade are oppo- 
site. Thus, considering one blade at a particular azimuth position ijr. 


= Ag cos p - (x - |)p - Jig sin p cos + 
^ X cos ilf + 2- X sin if 

a a 

and, for the opposite blade at ijr + it, 

= Ag cos P + (x - I)p + Jig sin p cos ^ + 

Q *n 

^ X cos •^ + g. X sin ijr 


Solution of Blade-Flapping Equations 

Blade-flapping an^Utude and stability may be studied by solving 
differential equation (19) or (20) for the flapping amplitude p as a 
fimctlon of azimuth angle ijr. Because the equations are nonhomogeneous 
and nonlinear with varying coefficients that contain integral expressions 
idilch became very con^lex when such -items as stall anfl con^pressibility 
are being considered, an esqpllcit solution is not feasible anfl numerical 
step-by-step methods must be used. 

The suggested procedure for solving equation (19) or (20) may be 
outlined as follows: 

(1) At some convenient azimuth position, assume arbitrary initial 
ve^Aies of p = p = 0. 

(2) Con5)ute the thrust moment at this selected azimuth position by 
evaluating the thrust-moment integral at several radial stations and 
integrating the results by numerical or graphical means. 

(3) Solve the resulting equation for p. 

(4) Bepeat the above steps at the next succeeding aziimxth position 
using v^ues of p and P obtained by integrating the quantity p. 

Q?he Bunge-Kutta method described in reference 6 is suitable for this 
phase of the problem. 
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-(^) Coirfalnue 1;lils process at siiccesslve azimuth, positions imtll 
P either reaches a steady-state value (that is, repeats itself periodi- 
cally) for a stable condition or diverges for an unstable one. 

Some saving in confuting time is possible, particularly -when steady- 
state values are being cosnputed, by using estl^ted values of p and p 
in step ( 1 ) rather than assuming them equal to zero. 

Physically, the blade response as calculated by this method may be 
visualized either as that resulting from the transition of the rotor 
blade from an arbitrary out-of-trim position to one of equilibrium (for 
a stable case) or, alternatively, the motion resulting from the sudden 
release of a blade that has been restrained from flapping. 

The method is also applicable to computing the transient response 
of blade, flapping due to arbitrary control displacements (by merely 
inserting any desired cyclic- or collective -pitch control sequence into 
eq. ( 16 ) at some time after steady-state conditions have been attained). 
In addition, the method may be \ised to determine steady-state flapping 
values at fll^t conditions outside the scope of conventional rotor 
theory. 


APPUCATICaJS 


In order to illustrate some applications of the method that are of 
current interest, the flapping stability of a nmber of unloaded-rotor 
configurations operating at extreme tip-speed ratios, as well as the 
transient flapping motion occurring during a hl^-speed pull-up maneuver 
of a conventional helicopter, was investigated. Most of the cases that 
were studied involved the free-to-cone rotor, but two seesaw-rotor 
exan^les are included for con^arlson in the unloaded-rotor application. 
Aside from type of rotor, the primary variables considered are the tip- 
speed ratio pg and the blade mass constant 7 ' . 


Calculation Procedure 

The blade-flapping equations were evaluated by the Rvmge-IQitta method. 
Azimuth intervals of 20° were \ised, except where noted. (Althovi^ the 
Runge-Kutta method entails caiculations at intezmediate stations and is 
thvis more lengthy than some other numericeil methods that were tried, it 
was found to be more accurate.) The thrust-moment contribution at six 
radial stations was evaluated in determining the thrust-moment Integrals 
in equations ( 19 ) and (20) . The calculations were performed at the 
Langley Laboratory by the Bell Telephone Laboratories X-66744 relay com- 
puter, and approximately ^0 minutes was required for the calculations at 
each Eizimuth position for the free-to-cone rotor. In terms of the more 
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widely used IBM Card Programmed Electronic Calcvilator (with, a floating 
decimal board), the estimated time per station would be reduced to 20 min- 
vrbes. With one of the hi^er speed, machines, such as the Remington Rand 
Universal. Automatic Computer (UMIVAC) or the Ml type 7OI electronic data 
processing machine, the machine time per station woxild probably be reduced 
to somewhat less than 1 minute. The seesaw-rotor calculations require 
about twice the time as those for the free-to-cone rotor. 

The section lift and drag cuirves used in the calculations are shown 
in figure 2 and are plotted for an angle -of -attack range from 0 ° to 360°. 
This method of handling the angles of attack is particularly convenient 
for aizbomatlc-machlne ccmq)utatlon in that the reversed-velocity region 
is handled in the same manner as the rest of the disk and thus the use of 
special instructions is avoided. In order to make the calculated section 
angle of attack positive over a 360° range, 360° should be added to o^, 

whenever oc^, as calculated by equation (16), is negative. Thus, for 

exanq)le, a negative angle of attack of - 30 °, such as mi ght occur in the 
reversed-velocity region, would, by adding 360°, be eacpressed as 330°. 
Values of cj and below the stall axe represented in figure 2 by 

the customary linear lift curve and three-term drag polar, whereas values 
above the stall are approximately in accord with wind-tunnel tests of air- 
foils throu^ 180° angle of attack (ref. 7 ) • Con^jresslbility effects were 
not considered in the present examples, althou^, as previously mentioned, 
the method of analysis is capable of handling such factors. In addition, 
the blades had zero hinge offset, were untwisted, and operated with zero 
cyclic pitch, except where noted. 


Flapping-Stability Ejcanples 

Free-to-cone .rotor .- The calculated blade flapping angles at a flight 
condition corresponding to 6ui unloaded rotor at (Xg = 1.0 f 0 „ = 0°j 

V \ • I ✓ 

TVg = 0.038) axe shown in figure 3 for three different mass factors. In 

each case the flapping motion is stable, in that steady-state values axe 
either reached or approached in a few revolutions, althou£^ some initial 
overshoot of the final, value is evident for the li^test blade (7’ = 2.62) . 
As an indication of the effect of blade mass factor on flapping, the maxi- 
mum steady-state aaplitudes shown in figure 3 axe plotted in figure ij- as a 
function of blade mass factor. The Tnaxl.mum positive flapping amplitude is 
shown to Increase almost linearly with 7'; whereas the TnaylTmim negative 
angle 'is, for this particular case, almost independent of 7' . 

The effect of blade mass constant on the flapping stability and 
amplitude of a rotor operating at a miore extreme fll^t condition 
(lie = 2 . 2 ; 0 ^ = 0 °; ?ve = 0.029) is shown in figure 5 * For ease of 

representation, the vertical scale is varied for the various values of 
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7’ in order to account for the larger flapping angplitudes of the lifter 
blades. The motion is seen to be stable for all the values of 7* repre- 
sented, althou^ the flapping changes from predominately first hemnonic at 
the lowest values of y* to a prominent one-half harmonic at the hipest 
values of 7* . The deviation of the dominant frequency of the transient 
flapping motion from the more normal first -harmonic motion is discussed 
in reference 2 and may be attributed to the varying spring and damping 
terms in the blade-flapping equation of motion. The curves in figure 5 
show that the heavier blades take a longer time to reach the steady-state 
condition than the lifter ones do and thus may be considered to be more 
sluggish in response to a disturbance. Figure ^ also shows that the 
steady-state flapping anplitudes of the lifter blades appear -feo be 
impractical ly large . 

In order to determine the effect of azimuth-interval spacing on the 
res\ilts, the calculations for three values of 7' in figure 5 were car- 
ried out in two parts, the first part with 40 ° increments and the second 
part with 20 ° increments. It appears that the larger interval would yield 
about the same resiilt as the smaller one insofar as the stability of the 
motion is concerned and that the smaller interval would probably yield a 
more acc\rrate result for the anplitude of the steady-state motion. 

Flapping stability at an even more extreme flight condition (|is = 3*0 
0 ^ = 2 °j = -0.0636) is investigated in figure 6. The motion is seen 

to be stable for 7' = 0.10 aaad 0 . 42 , althou^ for y* = 0.42 the motion 
includes a one-half -per-revolution harmonic inqposed on the basic one-per- 
revolution frequency. For a lighter set of blades (7 = 2 '. 62 ), the motion 
is purely divergent in less than 1 revolution. 

Seesaw rotor .- Coniputations were made for .two seesaw rotors involving 
different flight conditions and blade mass factors and are shown in flg- 
Tire 7 « As in the free-to-cone configuration, the heavier blades show a 
relatively slow response; whereas the lifter blades show the character- 
istic initial overshoot. It should be noted that the fll^t condition 
coiresponding to the light (7' = 2.62) seesaw blades in figure TC^) is 
identical to that of the light (7’ = 2.62) free-to-cone blades in fig- 
ure 6. Althou^ the flapping amplitude of the seesaw blades is large 
for this condition, it is significant that the motion is stable; \diereas 
the flapping motion of the free-to-cone blades in the Identical fli^t 
condition is highly divergent. It thus appears that the present method 
of analysis could be \ised to indicate likely avenues for improvement . 


Transient-Response Example 

The calculation of transient response of a rotor blade to arbitrary 
control inputs \inder conditions involving extreme rotor stall was carried 
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oirt by this method for the case of a free-to-cone rotor operating at a 
tip-speed ratio equal to 0.3 • The helicopter tos assumed to be under- 
going a pull-up maneuver •which involved t'wo arbitrary sequences of longi- 
tudinal cyclic- and collective-pitch control movements. For the flight 
condition chosen, the rotor is only partially stalled inboard before the 
start of the pull-up and undergoes extensive stall over most of the 
retreating side of the disk during the maneuver. 

The results of the calculation are plotted in figure 8, in which are 
shown both the blade-motion and control-movement time histories. The 
steady-state an^litude prior to the control motion ■was calculated as in 
the preceding cases. The steady state is very quickly reached (-within 
2 revolutions) and the amplitude is reasonable. When the cyclic- and 
collective-pitch controls are moved in accordance -with a pull-and-hold 
procedure, the blade motion responds almost immediately and reaches 
a new steady-state amplitude that is considerably larger tha n the 
anq)lltude existing before the pull-vq). If, instead of a pull-and-hold 
procedure, the controls are partially returned to their equilibrium posi- 
tion, the mavimiim positive flapping is reduced, as would be expected, 
with but little effect on the mpyimum negative amplitude. It •will be 
noted that the mavimum negati-ve anplitude is large and that the method 
of analysis permits calculations to be made for the study of blade- 
clearance problems in maneuvers. Similar calculations can be made to 
investigate the effects of varying the rate and amount of control motion, 
the phasing between the cyclic- and collective-pitch controls at ■various 
flight conditions, and the. effect of the pitching^ or rolling velocity 
developed diaring the maneuver. 


CQNCHJDING RMARKS 


A numerical method has been presented for studying the transient 
flapping behavior as well as for calculating the steady-state blade flap- 
ping motion of lifting rotors operating at extreme conditions involving 
stall and compressibility effects and at very hi^ tip-speed ratios eind 
rotor an^es of attack. The method can be applied to free-to-cone aufl see- 
saw rotors and to blades of any airfoil section, mass distribution, t-wist, 
plan-foim taper, root cutout, and flapping-hinge geometry. 

The method -was applied to the investigation of the effect of tip- 
speed ratio (equal to or greater than 1.0) and blade mass factor on the 
flapping stability of vinloaded free-to-cone and seesaw rotors the 
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effect of a pull-up maneuver on tlie 'transient blade response under con- 
ditions involving extreme stall at a tip-speed ratio equal to 0.5* 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 3^ 195^* 
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(b ) Velocities and angles at .blade element . 

Figure 1.- Forces and velocities acting at a blade element. 
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Figure 2.- Concluded. 
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